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ABSTRACT 

Weak gravitational lensing changes the angular power spectra of the cosmic microwave back- 
ground (CMB) temperature and polarization in a characteristic way containing valuable in- 
formation for cosmological parameter estimation and weak lensing reconstructions. So far, 
analytical expressions for the lensed CMB power spectra assume the probability density func- 
tion (PDF) of the lensing excursion angle to be Gaussian. However, coherent light deflection 
by nonlinear structures at low redshifts causes deviations from a pure Gaussian PDF. Work- 
ing in the flat-sky limit we develop a method for computing the lensed CMB power spectra 
which takes these non-Gaussian features into account. Our method does not assume any spe- 
cific PDF but uses instead an expansion of the characteristic function of the lensing excursion 
angle into its moments. Measuring these in the CMB lensing deflection field obtained from 
the Millennium Simulation we show that the change in the lensed power spectra is only at 
the 0.1%-0.4% level on very small scales (Afl < A' , I > 2500) and demonstrate that the 
assumption of a Gaussian lensing excursion angle PDF is well applicable. 

Key words: cosmology: large-scale structure, gravitational lensing, methods: analytical, nu- 
merical 



1 INTRODUCTION 



X 



Cosmological parameter analysis as well as investigation of inflationary models require a precise knowledge of the fluctuations of the cosmic 

microwave background. Therefore supreme technical eff'orts are made to equip current and future CMB experiments with instruments with 

. , ever greater resolution including sensitivity on CMB polarization. In order to interpret these low-noise data correctly a profound understand- 

frt ' ing of the physics of the primary CMB is not longer sufliicient but requires an equivally good understanding of CMB foregrounds. 

" " ' One of the most important of these is weak gravitational lensing: on their way from the last-scattering-surface to today's observer 

CMB photons are deflected by the intervening large-scale structure. Although each single deflection is small , their cumulative effect changes 
the statistics of the CMB fluctuations and accordingly their power spectra observed on today's sky (see iBartelmann & SchneideallOOll : 
iLewis & Challi noi^ 200ff). At prese nt, the le nsing effect has been detected in the temperature data of the WMAP satellite at moderate 
significance ([Smith etal. 2007; Hir ata et al.l 12 008). Apar t from being a CMB foreg round, th e lensing pattern itsel f is a valuable source 
of co smological information (Guziketal. 2000; Hu & Okamoto 20021; iHirata & Sell aki 120031 ; lOkamoto & Hij|2003l ; lAmblard et alj|2004 
iLessourgues et al 2006.;.Perotto et al._2006.;.de Putter et al...2009.). 

SeliakI <1996l) ; IZaldarriaga & SeliakI ( 119981) ; IChallinor & Lewia ( |2005[) developed efficient and accurate methods for computing the 



lensed CMB power spectra starting from the corresponding correlation functions. These methods all work under the assumption of a Gaussian 
distributed lensing excursion angle, which measures the difference between the deflection of two nearby light rays. Since the deflection angle 
is given by the gradient of the lensing potential one commonly accounts for the impact of nonli near structures on the lensed CMB spectra 
semi-analytically by applying nonlinear but Gaussian corrections from the HALOFIT model of ISmith et alj 1 120031) to the lensing potential 
power spectrum. This approach takes the additional power on small scales due to nonlinear evolution into account but ignores that coherent 
deflection by nonlinear structures at low redshifts alters the Gaussian character of the lensing excursion angle distribution. In this work we 
investigate how exactly this non-Gaussianity influences the lensed CMB power spectra. 

This paper is structured as follows. In Section|2]we develop the formalism for computing the lensed CMB temperature power spectrum 
in case of a non-Gaussian lensing excursion angle PDF. We relegate the discussion of the power spectra involving CMB polarization to 
Appendix lAl Section[3]is devoted to our numerical results. Here we quantify the non-Gaussianity present in the lensing excursion angle PDF 
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which we compute from the CMB deflection field of the Miflennium Simulation (MS). With this deflection field we then derive the lensed 
CMB power spectra and compare our results to power spectra resulting from using a pure Gaussian lensing excursion angle PDF. In Section|4] 
we summarize our results and give an outlook on future investigations. Finally, in Appendix |B] we give explicit expressions for the lensed 
CMB correlation functions used in our numerical computations. 

Throughout this work we assume a spatially flat ACDM cosmology with adiabatic Gaussian initial perturbations. The relevant parameter 
values are: n„, = 0.25, Ha = 0.75, Hq = 73kms"' Mpc"', Q.i, = 0.045, n, = 1.0, ag = 0.9 and r = 0.0 (no primordial gravitational waves 
present). These parameters are equal to those of the Millennium Simulation JSpringel et al.ll2006h . 

2 FORMALISM 

2.1 Lensed CMB temperature power spectrum 

Working in the flat-sky limit the 2D lensed temperature field is given by the remapping 

0(x) = 0(x + a(x)) (1) 

mediated by the lensing deflection angle a, i.e. the gradient of the lensing potential if/: ex = Vifi. In the approximation of instantaneous 
recombination the CMB can be described by a single source plane at conformal distance x = ;t ■ In the absence of anisotropic stress the 
lensing potential is then given by the line of sight projection of the physical peculiar gravitational potential (p: 

x' 

>p{n) = \ fdx^^-^<p(xn;x), (2) 



where the underlying geometry is flat JBartelmann & SchneideiiuOOlc iLewis & Challinoiil200q) . Introducing the Fourier transform of the 
temperature field via 

0(x) = j ^0(iy'", (3) 

the spectrum for a statistically homogeneous and isotropic field reads 

{0(i)0*(i')) = cf%(\ - n (4) 

Then, ignoring the weak large-scale correlation between CMB and lensing potential due to the integrated Sachs-Wolfe effect, the lensed 
correlation function of the CMB temperature fluctuations is given by 

|(r) = (0(x)0(x')) = r -^Cf^e"" (e'K^W-^c')]) (5) 

^ ' J (2;r)- ^ ' 

where r = |x - x'|. It is worth noting that the lensed correlation function only depends on the relative displacement, the so-called lensing 
excursion angle, 5a(r) = q;(x) - q;(x') and that this dependence is given by the characteristic function of the lensing excursion angle: 



/- 



<^5c.(1) = (e""") = J d(<5a) p(Sa.)e'''"- . (6) 

The last equality reveals that the characteristic function is the Fourier transform of the PDF p (Sol). Hence, it carries the same information as 
the PDF itself. From equation lO the lensed power spectrum is readily obtained by 

Cf'^ = 2n ^drrl(r)Mlr) (7) 

where J„{z) denotes the n-th order Bessel function dAbramowitz & Stegunll 19651) . 

2.2 Lensing excursion angle 

In linear theory the lensing potential is a Gaussian field and so is its gradient, the lensing deflection angle. Accordingly, the lensing excursion 
angle is a Gaussian variate and therefore the characteristic function of the lensing excursion angle is just given in terms of the variance 



(8) 



^^u.[«(x)-«(x')]^ = exp|-i ([1 • SaA = exp(-i/2 [aHr) + cos 2(0, - cp,)C,,2(r) 

where ^i,. denotes the angle between 1, r and the x-axis. We have defined cr^(r) = USa^j = j (Cgi(O) - Cgi(r)j. Cgi and Cgi_2 are given in 
terms of the power spectrum of the lensing potential Cf JChallinor & LewisluOOSh : 

CAr)= — fdlfcf^MIr) and Cnir) = — f dU^Cf^Jiilr). (9) 

"" 2n J °' Itx J 

Insert ing equation ^ into equation (O and performing a perturbative exp ansion in C^\2 up to se cond order one recovers the expressions 
derived bv lChallinor & LewisI |2005). However. iHamana & Mellieg ( 1200 ih and lHamana et alj J2005r) showed in numerical weak lensing ray- 
tracing experiments that the lensing excursion angle is not Gaussian distributed. Its PDF has indeed a Gaussian core but also exponential 
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wings. These wings are a consequence of coherent scattering by individual massive haloes with mass larger than lO''' Mq//i. Since coherent 
deflection demands a sufficiently small intrinsic separation of the light rays the exponential wings appear prominent in the excursion angle 
PDFs obt ained from light rays with intrinsic separation of a few arcminutes and are negligible in those of separations larger than one degree 
JHamana et al. 2005). The contributions from coherent scattering broaden the PDFs, i.e. larger excursion angles are more probable than in 
case of a pure Gaussian PDF. For example, for intrinsic separations smaller than two arcminutes an excursion angle of one arcminute is about 
ten times more probable than for a Gaussian PDF. 

2.3 Non-Gaussian probability density function 

In order to investigate how the non-Gaussian features of the lensing excursion angle PDF affects the lensed CMB temperature spectrum one 
has to use (in principle) all moments of the lensing excursion angle for computing the lensed correlation function (|5}, which is now denoted 
with a hat to distinguish it from the correlation function derived under the Gaussian assumption (denoted with a tilde) 



(10) 



(11) 



For the last equality we used the binomial law 

(1 • Sc.)" = {l.Sa. + IM,) = g (J(/,.foO [IM,) with (J = ^^^— ^. 

Following iKamionkowski et al.l ( 119971) one should choose the local coordinate system, in which to define the correlation function, aligned 
with the great circle coimecting the two points where the temperature fluctuations are measured. In the flat-sky limit this choice of coordinate 
system corresponds to evaluating the two-point correlator at the origin of the flat coordinate system and at a point on the x-axis at distance 
r. One of the great advantages of this coordinate system is that here the correlation tensor of the lensing excursion angle is diagonal, i.e 
(aj{x)aj(x')\ oc Sjj. In case of a Gaussian distribution it then follows immediately by virtue of Wick's theorem that different moments of 
different components are statistically independent. In Section [373l we will show that it is reasonable to assign this property also to the (non- 
Gaussian) PDF of the lensing excursion angle. Hence, introducing polar coordinates 1 = (Zcos0, /sin0), the lensed correlation function 



reads 

i(r) = 



1, 1 1 d^ m cr^-' t'ht 0^" -^* ^ -""" ^ K-) {'<■') ■ (12) 



A further simplification can be obtained by demanding that components of the lensing excursion angle are distributed symmetrically about 
zero. This assumption is very natural since otherwise there would be a preferred direction along the coordinate axes. For a symmetric PDF 
all odd moments vanish, hence, 






„=0 k=0 q=0 

(13) 



To get the second line we used the fact that 

d0 e'-~'""1' cos" (j, = 2ff(-i)" -^ Uz) (14) 

dz" 

and that 

^ ' " tiz). (15) 



/ 






The expression for the lensed CMB temperature correlation function given in equation ( 1131 ) is exact for any PDF of the lensing excursion 
angle that is symmetric about zero and whose moments of different components are uncorrelated, which we have shown to be valid in the 
deflection field obtained from the Millennium Simulation (cf. Section [J!2t . 

3 NUMERICS 

3.1 Truncation 

For explicitly computing the lensed CMB temperature power spectrum via equation ( I13t one has to truncate the expansion of the characteristic 
function at a certain order n. To find a reasonable value for n we resorted to the approximation of a purely G aussian distributed l ensing 
excursion angle and verified the performance of the series expansion in comparison to the lensing method of IChallinor & Lewis! ( 120050 
described in Section |2J2] and numerically implemented in CAM^J. Successively increasing the order n taken into account in the series 

' http://camb.info/ 
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expansion we determined tliat n for wliicii botli lensing methods worlc equivally well. We confirmed that in the case of a pure Gaussian 
lensing excursion angle PDF for n = 3 almost perfect agreement between both lensing methods can be achieved. The deviations are largest 
on very small scales but do not exceed (9(10"''). Thus, for the numerical implementation of our lensing method we truncated the series 
expansion in equation l ll3t at n = 3, i.e. we included the sixth moments in the series expansion. Explicit expressions for n = 3 are given in 
Appendix IB] 

3.2 Moments 

Aiming at the influence of the non-Gaussian features in the lensing excursion angle PDF on the lensed CMB power spectra we ca nnot 
compute the moments needed for its calculation analytically using the power spectrum approach developed by ISeliakI ( 1 19941 . 1 19961) (cf . 



Section |X2t . This approach is well suited for describing gravitational light deflection by the intervening large scale structure but since it 
is based on linear perturbation theory it does not account for coherent lensing scatter by individual massive haloes, which gives rise to the 
exponential wings in the lensing excursion angle P DF (cf. Section | 2 . 2[>. The refore, we either have to resort to numerical simulations or to use 
the halo model of large scale structure reviewed bvi Coorav & ShethI j2002 h. 

In this wo rk the moment s neede d for the computation of the lensed CMB power spec tra are obtained from the lensing deflection field 
constructed bv ICarbone et al.l <2008l) which is based on the Millennium Simulation (MS) JSpringel et al.ll200q) . Being an all-sky map the 
deflection field is given as angular gradient of the lensing potential 

a(n)= Cfl— -l-e^^— — U(ii) = ag{n)e^e + a^{Vi)t^. (16) 

\ aO sm(pd(l>j 

Since the non-Gaussian features of the lensing excursion angle PDF are only substantial for intrinsic light ray separations smaller than one 
degree (cf. Section [Z2t . i.e on scales where the curvature of the sky is negligible, we approximate 

Sa(r) = a(x) - a(x') X a(n) - a(h') = Sa(fi) with r x fj and n ■ n' = cos/J, (17) 

where the bars indicate the basis defined by the geodesic connecting n and n'. The error of this approximation is slightly increasing with the 
light ray separation and finally reaches ~ 2% for two rays intrinsically separated by one degree. 

3.3 Correlation coefficients 

The expression for the lensed correlation function il3t assumes that different moments of different components of the excursion angle are 
statistically independent. To show that this is indeed the case we compute the correlation coefficient defined by 

{(X - {X}){¥ - {Y})} 

^J{ix-{x)r-}^{(Y-{Y)r-y 



P(X, Y) = -^ _] '/; ^^.^i^i— - , X,Y e{Sa';„i,\n = l,...,5\ (18) 



for all relevant combinations of X and Y. All p(X, Y) are compatible with zero at the level of 10 ', justifying the assumption of statistical 
independence. 

3.4 Non-Gaussianity 

Equation ( 113b reveals that all information about the non-Gaussianity of the lensing excursion angle PDF is carried by its moments. It is 
therefore natural to quantify the amount of non-Gaussianity of the excursion angle PDF via comparing its moments with those of a Gaussian 
PDF. The moments of the latter can all be expressed in terms of the variance: 

(^X^") = (2n - ly.l (^X^y with (2n - 1)!! = (2« - 1) ■ (2« - 2) • ... • 5 • 3 • 1. (19) 

Thus 

and 7/ = -^ (20) 



{r-y (r-y 

or rather their deviation from 3 and 15 carry information about the amount of non-Gaussianity present in the fourth and sixth moment of 
any symmetric PDF. Consequently, in order to ensure to capture only the effects on the lensed power spectrum which arise from the non- 
Gaussianity of the lensing excursion angle PDF it is recommended not to use directly the moments of the excursion angle components but to 
use instead /cg^ and T/e,^ defined by 

Ke.m^ ; ■'* ' and iie.^(J3) ^ j ^ (21) 

in analogy to equation MQ\ . Thus we used 

{Sal{r)) = a\r) + C^y2{r), {5a]{r)) = K„{r){Sal{r)f and {Sa'^(r)) = rio(r) {Sal(r)f (22) 

and 
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Figure 1. Quantifying the non-Gaussianity of the PDFs of the lensing excursion angle components, obtained from the MS deflection field, by the deviation 
of the fouilh (left panel) and sixth (right panel) moment, respectively, from the Gaussian case, i.e. Kg,/, = {Sa^ .) / [Sal A and rie,p = ( 5(5* , ) / {Sag , ) . The 
deviations from the Gaussian expectation (Kg,/, = 3, rjg^,/, = 15) decrease rapidly with increasing intrinsic light ray separation j8. For separations larger than one 
degree there are only small amounts of non-Gaussianity left. 



{Sa;{r)) = a\r) - Qu(r), (fo^(r)) = K^r) {Sal{r)f 



and 



{dat{r)) = u4r){s4('-)f 



(23) 



respectively, for our actual computation of the lensed correlation function accounting for th e excess of small sca le power of the lensing 
potential due to nonlinear growth via semianalytical corrections from the HALOFIT model (cf. lCarbone et alj2008h . For computing Kg^^p and 
J7e,(A from the CMB deflection field of the MS we sampled pairs of pixels with fixed angular separation and computed their differences in the 
local coordinate system defined by the connecting geodesic. The rotation angle needed for the transion from the coordinate system used by 
the MS (cf. equation[T6t to the geodesic basis can be readily obtained from identities of spherical triangles. From these samples we estimated 
the first three even moments of the lensing excursion angle components and computed /ce,,^ and rjg^^p. They are shown in Figure [T] Note that 
the sampling errors are negligible and are not shown. The non-Gaussian features are more pronounced in the distribution of the 0-component 
of the lensing excursion angle. The deviations from the Gaussian expectation are larger in case of t/b ,0 than for /fo^ as expected since the 
PDFs are broadened due the coherent scattering by nonlinear structures (cf. Section |Z21 (. Furthermore Figure [T] reveals that there is still a 
small amount of non-Gaussianity for light ray separations larger than one degree. The Gaussian expectation is not reached until fj ~ 20°. 
For such large light ray separations, however, the approximation Sotifj) a 8oL{r) given in equation l ll7t is not longer valid. Since the error 
from extending this approximation up to separations larger than one degree would fairly exceed the one resulting from neglecting the small 
amount of non-Gaussianity present for/? > 1° the latter is omitted in the remainder of this work. 



3.5 Results 

As seen in Section [X2l for intrinsic light ray separations up to one degree the lensing excursion angle can be well approximated by the 
diff^erence of two deflection angles defined in the geodesic basis. In the regime of larger separations the non-Gaussianity of the lensing 
excursion angle PDFs is weak and we confirmed that in this regime the Gaussian approximation is very good. Thus, it is natural to establish 
the following computation scheme for the lensed correlation function: For intrinsic light ray separations less than one degree one uses 
formula ( 113b together with the moments computed via equations ( I22t and J23t . For larger separations one uses equation ^ together with 
equation ([8j, which is valid for Gaussian distributed excursion angles. The same computation scheme can be applied to the lensed CMB 
polarization power spectra by using the corresponding equations given in the appendix. 

The influence of a non-Gaussian PDF on the lensed CMB power spectra can now be quantified by computing two sets of power spectra 
Cf^ and Cfj^f. For both sets the computation scheme for the lensed correlation functions described above is used but in case of the reference 
spectra we use Gaussian moments on all scales, i.e. setting Kejp{r) = 3 and r]g,p(r) = 15. Figure |2] shows both the lensed and reference CMB 
power spectra. In the lower half of each panel the ratio between the corresponding spectra is depicted. 

These ratios show that the influence of the non-Gaussianity in the lensing excursion angle distribution function is marginal. The dif- 
ferences caused by the exponential wings of the PDF are just several per mile on very small scales (A0 < 4', / > 2500). On large and 
intermediate scales all lensed power spectra but the 6-modes are almost unaffected. The lensed S-modes, however, are altered on all scales, 
reflecting the fact that in the cosmological model assumed in this work primordial gravitational waves are absent and thus polarization of the 
S-type is completely lensing induced. 

Taking into account that neglecting the curvature of the sky, as we did in this work, already causes an error in the lensed CMB power 
spectra at the 0.3% -1.0% level (Challinor & Lewis 2005) we conclude that the influence of the non-Gaussian part of the lensing excursion 
angle PDF on the lensed CMB spectra is far from being observable and can be safely neglected. Furthermore, it is very likely that uncertainties 
in the recombin ation history and contamination from various late-time secondaries affect the lensed power spectra in an even stronger way 
JLewis & Chall inor 2006). 

Our results have been anticipated by the purely numerical work of ICarbone et alj J2009l) . Using the MS they constructed lensed CMB 
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Figure 2. Lensed CMB power spectra computed via the method desciibed in the text accounting for the non-Gaussian features of the lensing excursion angle 
PDF discussed in Section lT2l (blue dashed curves). The reference power spectra (solid red curves) were computed by the same method but assume a Gaussian 
PDF of the lensing excursion angle. The ratio of the corresponding spectra is shown in the lower part of each panel (solid green lines). 



temperature and polarization maps and computed the corresponding power spectra. They agree well with the o nes obtained from CAMB 
including nonlinear corrections to the lensing potential from the HALOFIT mode l. Henc e, Carbone et alj ( 120090 increased the variance of 
the lensing excursion angle but did not drop the assumption of a Gaussian PDF (cf. lLewisll2005l) . 

Our analytical approach, however, takes the non-Gaussianity of the lensing excursion angle PDF explicitly into account revealing that its 
influence on the lensed CMB spectra is weak and therefore the assumption of a pure Gaussian lensing excursion angle PDF is well applicable. 
The weakness of the impact of the non-Gaussian features is due to the fact that they are only prominent on small scales (see Figure[T), while 
lensing smooths the CMB signal over a large angular range. 



4 SUMMARY 

The topic of this paper is a derivation of the fluctuation statistics of the lensed CMB temperature and polarisation, taking non-Gaussian 
features of the lensing deflection angle distribution into account. 

(i) Starting point of including non-Gaussian distributions of the deflection angle is the expansion of the characteristic function, which 
enters the computation of the lensed CMB spectra, into a a series in terms of its moments. For computing lensed CMB spectra with non- 
Gaussian deflection angle fields we provide a set of analytical expressions, and expand the expressions up to the sixth-order moment of the 
deflection angle distribution. 

(ii) Simulated deflection maps (provided by ICarbone et alj2009l) show considerable amounts of non-Gaussianity on small scales below a 
degree, which we quantified with the fourth and sixth moment. On angular scales of an arcminute, they exceed the Gaussian expectation by 
a factor of 1.5 and 3, respectively, and drop close to their fiducial, Gaussian values close to one degree. 

(iii) We show that deviations in the angular spectra relative to those derived assuming Gaussian statistics are most important on small 
angular scales, but remain below the percent level, confirming that the Gaussian approximation is very good. In particular we confirm that 
errors introduced into the spectra by non-Gaussian deflection angle statistics are smaller than those caused by other secondary anisotropies. 

(iv) One can give a simple physical argument, why the impact of the non-Gaussian lensing excursion angle PDF is so weak: On the 
subdegree scale, where the non-Gaussianity is considerabl e, the CMB spectra are almost featureless. As lensing cannot generate features in 
a featureless CMB (see lHiJuOOduLewis & Challinonl2006h the influence of the non-Gaussianity of the lensing excursion angle, which itself 
is only substantial on subdegree scales, is very weak. 
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(v) Estimates of cosmological parameters and weak lensing reconstructions are not seriously impeded by non-Gaussianities in the deflec- 
tion angle distribution, given the small differences relative to spectra derived with a Gaussian approximation. Furthermore, the difl"erences in 
the spectra involving the temperature and E-type polarization are substantial only at high multipole order where the signal-to-noise ratio is 
small and where the spectra do not possess strong parameter constraining power due to Silk-damping. 

The formalism presented here can be applied to investigating primordial non-Gaussianities by CMB lensing e.g. in terms of the /iML-model 
or the ^^-model, for which the higher-order moments are directly calculable. It would be interesting to see if lensed CMB spectra are 
significantly distorted in these cosmological models, because they inherently are able to provided stronger non-Gaussian features on larger 
angular scales compared to non-Gaussianities generated by nonlinear structure formation. 
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APPENDIX A: POLARIZATION 

A I Correlation functions and power spectra 

The polarization of the CMB is described by the Stokes parameters Q and U constituting the spin-2 polarization field P = Q + lU, whose 
decomposition in gradient-like E- and curl-like S-modes reads (in the limit of a flat sky) 

P(x) = - I — (£(1) - iB(l))e-"V''' (Al) 



■/s< 



JChallinor & Lewisll2005h . Defining the correlation functions involving P in the same local coordinate system as before, i.e. with the x-axis 
adapted to the vector connecting x and x', we have 
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Ur) = {P'{x)P{x')), Ur) = {P{x)P{x')) and ^,(r) = <0(x)/'(x')> 
JChallinor & Lewia2005l) . The corresponding power spectra are then calculated via 

Cf^ + Cf^ = 2n \ rdr Jo{lr)^+{r), Cf^ - Cf = In \ rdr Ji{lr)^-{r) and Cf^ = 2n f rdr JiW^xir) 

JLewis & Challinoill2006h . 



(A2) 



(A3) 



A2 Lensed correlation functions 

The lensed correlation functions involving the polarization field valid for a general distribution function of the lensing excursion angle can 
be derived in complete analogy to Section l23] The expression for i+(r) is identical with equation ( 113b if one replaces Cf^ by Cf^ + Cf^: 

(il ■ da)" 






n-k l(k+ci) s„+k+r , 



2^ J ^ '^riU^rihi (2n)!4*+H 



(-1 )"+*+>■ /2n\/n - k\l2{k + q) 



i,=0 k=0 i:/=0 ;-0 



2k l\ q 



J2ik.,-r,(lr){Saf){6a^!-"-'>). 



(A4) 



Some additional effort, however, has to be put in the computation of i-(r) and ix{r), since here the additional factors of e ^"^ in the spin 2 
polarization do not cancel. After a somewhat lengthy but straightforward calculation we find 



^-''' = /(S^(^'^-^'^)-^4^^"i|- 



(il ■ Sa)" 



11=0 



I2(k + q+l)' 



±/M,(cr-C,")gggi^g)(-')(..f)K-») 

2(k+q+2) , , 2(t+<;+l) , 

■ 8 Yj ^-^A ^ l/2(t+,+2-.)(/'-) - 8 Yl (-l>t 

r=0 \ '' I r=0 \ 

2{k+q) I 



J2(k+q+l-r){lr) 



J2(k + q) 



and 

Ur) 



J2(k+q-r)(.lf) 



(il • 8a.)" 






11=0 k=Q q=Q 

2{k+q+i) 



!4*+«\2yt/\ g 



2 2 (-1)' 



/2(^ + ^+l) 



2(*+f/) 



J2ik+q+\-r){lr) + ^(-1)' 



I2ik + q) 



J2{k+q-i-)(.lr) 



(A5) 



(A6) 



(A7) 



APPENDIX B: EXPLICIT EXPRESSIONS 

Here we give explicit expressions for the lensed correlation functions used for our numerical calculations. Truncating at « = 3 in the 
equations l |13l l, l lA4t . l |A5l l and l |A7t , we obtain: 



l{r) = ^ J/d/ Cf^ijodr) - U^ [((fcj) + (6aff) Mir) - ((fo;) - (Saj)) Jiilr)] 
^ {{6at) + {Saf}) (3Mlr) + Mir)) - i ((fo^) - (fo^)) Mir) 



+ — f 
24 



720 



+ -{6af){6a]){JSr) - JSr)) 
-^ ((fc?) + {dafj) (5MIr) + 3MIr)) + ^ {{Sa^) - (Sa^)) (-l5MIr) - Mir)) 



15 



— (Sal) {6a*) {2Mlr) + Mir) - 2Mlr) - Mir)) 



32 



15 

'32 



[da*,] (dal) (2Mlr) - Mir) - 2J^{lr) + J(,{lr)) 



(Bl) 



f+('-) 
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^ {{6al) + {6a^)) (3Mlr) + Mir)) - )- {{Sa'*^) - (Sa^)) J2(lr) 



+ — f 
24 



720 



+ -{6a'^{5al){Ulr)-J^{lr)) 
-^ ((fc^) + (6al)) (5Jo(lr) + 3Mlr)) + ^ ((fc«) - (fo«)) (-IS/jC/r) - Mir)) 



15 

'32 
15 

^32 



'6al) (Sa^) (2Mlr) + JiUr) - 2J,{lr) - Mir)) 
(Sa'^) {6afl UMlr) - Mir) - 2J,ilr) + Mir)) 



+ ... 



Ur) 



+ — f 
24 



i^ JzdZ [Cf - Cf ) [jSr) - \f {{Sal) + (dal)) Mir) - 1 ((fcj) - (Saj)) (Mir) + Mir)) 
^ ((fo^) + (Sa^)) (Mir) + 6Mlr) + Mir)) - j ((fcj) - (fo^)) (Mir) + Mir)) 
-^ (Sal) [Sa]) (Mir) - 2],(lr) + Mir)) 



720 



-7T {{Saf} - (^o'v)) a^Mlr) + ISMlr) + Jiodr)) 



11 
"64 

15 

^64 



'6(4) (Sa^) (-2Mlr) + AJ,(lr) + MM - 2Mlr) - Jw(lr)) 
{Sa'^) {Safl (-2Mlr) + 4Mlr) - Mir) - 2Mlr) + Jio(lr)) 



and 

ix(r) 



i- Jzd/ Cf^Mlr) - ^f {{6al) + (Sal)) Mir) - ^ ((fo^) - {Saf}) (Mir) + Mir)) 



1 

/ 

24 



+ ^1" 
720 



^ {{Sat) + {6af}) a Mir) + Mir)) + j {(Saf) - {&af)) (Mir) + Mir)) 



16 



-{6al){6aj)(-Mlr) + Mlr)) 
^ ((fo«) + {6a^) (ISMlr) + 3Mlr)) 



32 



+ 7T {{Saf} - (Sc^y)) aSMlr) + 16Mlr) + Mir)) 



11 
"64 



[6a^ (Sa^) (-Mir) - 2Mlr) + 2Mlr) + Mir)) 



15 



— (Sat) (daj) (Mir) - 2Mlr) + 2Jdlr) - Mir)) 



64 



(B2) 



(B3) 



(B4) 
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